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We present the analytical calculation of entanglement entropy for a class of two dimensional 
field theories governed by the symmetries of the Galilean conformal algebra, thus providing a rare 
example of such an exact computation. These field theories are the putative holographic duals to 
theories of gravity in three-dimensional asymptotically flat spacetimes. We provide a check of our 
Held theory answers by an analysis of geodesics. We also exploit the Chern-Simons formulation of 
three-dimensional gravity and adapt recent proposals of calculating entanglement entropy by Wilson 
lines in this context to find an independent confirmation of our results from holography. 
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Introduction. Entanglement — originally coined 
“Verschrankung” by Schrodinger [1] - is one of the most 

fascinating and mysterious features of quantum mechan- 
ics. It has been the basis of development of new branches 
of physics like quantum information and communication 
[2]. For these applications, as well as for theoretical anal- 
yses, it is useful to quantify the amount of entanglement 
of a quantum system, and various measures have been 
proposed, see e.g. [3]. From a theoretical perspective par- 
ticularly the entropy of entanglement, or entanglement 
entropy (EE), has emerged as a valuable tool. 

The calculation of EE for interacting quantum field 
theories (QFTs) remains a daunting task. One of the few 
examples where analytical calculations can be performed 
is the EE in two-dimensional (2d) conformal field theories 
(CFTs) [4]. Here one can invoke the infinite dimensional 
symmetries of the Virasoro algebra — the algebra that 
generates conformal transformations in 2d — to simplify 
the calculations. E.g. the EE of a strip of length £ of a 
2d CFT with equal left and right central charges c on the 
infinite line is given by [4-6] 

^ E ° E FT2 =|ln- (1) 

3 a 

where a is the lattice size that regulates the ultraviolet 
divergences. The problem of calculating EE in QFTs 
quickly becomes intractable with increasing dimension 
and complexity, since in general there is no infinite di- 
mensional symmetry algebra to fall back on. 

The holographic principle has played a vital role in the 
recent resurgence of EE calculations for QFTs. Accord- 
ing to this principle, a theory of (quantum) gravity in 
a certain spacetime is completely equivalent to a QFT 
without gravity living on the boundary of that space- 
time. The celebrated Anti-de Sitter/conformal field the- 
ory (AdS/CFT) correspondence has given this bold pro- 
posal a firm footing in the context of string theory by re- 


lating Type IIB superstring theory on AdSsCgiS 5 to A/ - = 4 
supersymmetric Yang-Mills theory on its boundary [7]. 

As a remarkable consequence of holography, the calcu- 
lation of EE in CFTs has been conjectured to be equiv- 
alent to the computation of the area of an extremal co- 
dimension-two-surface in AdS [8] . There are recent devel- 
opments which have built towards a proof of this relation 
[9-11]. Thus, using holography, the involved calculation 
of EE in higher dimensional CFTs has boiled down to the 
fairly simple computation of areas of extremal surfaces in 
AdS. See [12-27] for selected pre-cursors and elaborations 
on the connection between gravity and entanglement. 

The present work has two main goals: to calculate EE 
for a novel class of QFTs, and to derive EE for the same 
class of QFTs holographically. 

Related to the first goal, we shall calculate EE for 
Galilean conformal field theories (GCFTs), thus provid- 
ing a rare example of an exact computation of EE. 

Concerning the second goal, 2d GCFTs [28] were pro- 
posed as duals [29] to three-dimensional (3d) gravity in 
asymptotically flat spacetimes [30] . There has been con- 
siderable interest in trying to understand holography in 
flat spacetimes in the light of these findings [31-49]. We 
shall provide further evidence for a flat space/ G CFT 
correspondence by holographically calculating EE in 3d 
asymptotically flat spacetimes, thereby confirming the 2d 
GCFT calculations. 

EE in 2d CFTs. Let us now briefly recall the strategy 
employed for the calculation of EE in 2d CFTs, starting 
with the definition of EE. If p is the density matrix of a 
quantum system with multiple degrees of freedom, and 
the Hilbert space of the system can be written as the 
direct product TL = TLa®'Hb, the reduced density matrix 
Pa for the subsystem A is defined by the partial trace over 
the Hilbert subspace TLb, Pa = Tr Bp- The EE is then 
the corresponding von Neumann entropy 

S A = -Tr (p A In pa) • 


(2) 
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EE is conveniently calculated by employing the replica 
trick, using Renyi entropies, S A ^ = y^lnTrp^. EE (2) 

emerges as the limit Sa = lim„_n S A . More specifically, 
in a 2d QFT defined on a lattice, where the subsystem 
A consists of n disjoint intervals, one computes Renyi 
entropy S A ' by making n copies of the system with one 
disjoint interval and sewing them together. This defines 
an n-sheeted Riemann surface with a partition function 
Z n (A) and Trp^ = z"(A) n • From this one can calculate 

the EE of the system as a limit Sa = — lim„^! ^ Zi(A) n • 
For details of this construction see e.g. [50, 51]. 

We focus here on the EE of a single interval [u, v\ of 
length £ = |u — u| in an infinitely long one-dimensional 
quantum system at zero temperature and recapitulate 
the main insights and results. One key ingredient is to 
map the n-sheeted Riemann surface 1Z„, to the complex 
plane C, where n labels the n th Renyi entropy. An- 
other relevant aspect is that in a 2d CFT, the energy- 
momentum (EM) tensor on TZ n is related to that on C 
by T(w) = (Jj|) 2 T(z) + j^{z,w} where c is the central 
charge of the CFT and {,} the Schwarzian derivative. 
Taking expectation values of both sides and exploiting 
(T(z)) = 0 on C yields 

Wk. = (l - W ( V >a W < 3 > 

24 V n z / (w — u) z (w — v) z 
Conformal Ward identities then imply 


(T(w)$ n (u)®-n(v))c 

($ n (u)$- n (v))c 


( 4 ) 


where are twist fields in the 2d CFT, viz., pri- 

maries with conformal weights A = A = (c/24)(l — 
1/n 2 ). The key observation is that Trp^ transforms un- 
der a general conformal transformation as a 2-point func- 
tion of such primaries. From this, Calabrese and Cardy 
deduced Tr p A oc c( " 1 ' / ' 1 ^ 6 . The ensuing expres- 

sion for Renyi entropy 


Q (n) 
D A 



( 5 ) 


[M n , M m ] =0 n,meZ ci,c M el. (6) 


We consider now the representation theory of this 
Galilean conformal algebra (GCA) to construct the 
Hilbert space of a GCFT. We choose to work with the 
highest weight representation [52, 53]. Any state in 
this field theory is labeled by two weights = 

hL\hL,h M ), M 0 \h L ,h M ) = h M \hL,h M )- Following CFT 
literature, we define primary states by demanding that 
they be annihilated by all generators L n , M n with posi- 
tive n. We shall use the plane representation of the GCA 
(6) 

L n = x n+1 d x + (n + 1 )x n y d y , M n = x n+1 d y (7) 

where x , y are our 2d coordinates. We define the EM 
tensor of a GCFT on a plane [54] by 

T A) = Z( i " + ( n + 2)f M n )x- n ~ 2 (8) 

n 

T(A = y] M n x~ n ~ 2 . (9) 


An important result is that the 2-point correlator of pri- 
mary fields {$ l (xi, yi)} with weights h l L ,h l M has a par- 
ticularly nice form in the plane representation, fixed by 
the symmetries of the GCA (up to normalization): 

(® 1 {x 1 ,y 1 )® 2 (x 2 ,y 2 )) ocx]" 2 2hi e~ 2hM ^ 6^ 5 h i (10) 

n L rL M 

Here and below expressions like X \2 always denote differ- 
ences between coordinates, X 12 = X\ — X 2 - 

The 3-point correlator of the EM tensor with two 
primary fields (of weight {/i£,/iA/}) is determined from 
Ward identities and the 2-point correlators (10): 


(' T {1 ){x w ,y w )$(xi 1 y 1 )<S>(x 2 ,y 2 )) = ) (11) 


/ y 12 _ Vw\ _ Vw 2 \ 
\Xi2 X w i X w 2 ' - 


Xw\Xyj2 ' 

—2h L -2 h M \ 
l 12 e 


tiL — 2 h M \ 

(T(2)(xw,yw)®(x 1 ,y 1 )$(x2,y 2 )) = (—^—) ( 12 ) 


7 — 2hL — 2hjvi 

h M x 12 e “is 


in the limit n — > 1 then establishes the result (1) for EE. 

GCFT symmetries and correlators. GCFTs are the 
non-relativistic analogues of CFTs. They can also be 
defined independently from any relativistic parent CFT 
[28] . We now build towards our goal of calculating EE in 
GCFTs. A first step is to analyze the symmetry algebra, 
its highest weight representations and restrictions from 
Ward identities on correlation functions. The symmetry 
generators L n , M n of a 2d GCFT obey 

[An? L m ] — (n ??z)Tnym T — [jl 7z) (in+m, 0 
[_L n , M m ] — (ji m)M n+m + ^ (ji n) ^n+m, o 


The 3-point correlators above will play a crucial role in 
our derivation of EE for GCFTs. 

EE in 2d GCFTs. We are now ready to construct the 
EE for a GCFT of an infinite one-dimensional system for 
one interval at vanishing temperature. 

We intend to use the analogue of (4) in 2d GCFTs 
in order to be able to make the statement that also in 
GCFTs Tr p A depends only on the form of the 2-point 
function of primary operators. We have already com- 
puted the right hand side of (4) from symmetry consid- 
erations above. For arriving at the left hand side, we 
shall employ a limiting technique from 2d CFTs [71]. 
We need the GCFT equivalent of (3), assuming a smooth 
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non-relativistic limit. The expressions for the GCFT EM 
tensor in terms of the original CFT EM tensor read: 
T(i)(x w ,y w ) = lim e _> 0 [T(w) + T(w)] , T (2) {x w ,y w ) = 
lim e _,.o e [T(w) — T(w)] with w = x w + ey w and the 
limit is defined by e — »• 0. Denoting u = X\ + ey\ and 
v = x 2 + ey 2 , this leads to: 


G(i) (x w , Vw))n n 


G(2) (x w , Uw^TZrt 



1 \( 

' Xi 2 



n 2 )\ 

- Xyj\X w 2 

) 

C L 

cm 

ryvi 

Uwl 

24 

12 

VcCl2 

Xw 1 


J_\ /_^12_\ 2 £m 
n 2 J\x w i x w2 J 24 



(13) 
fc2V 

(14) 


Combining the results (10)-(14) obtains (i = 1,2) 


(T(i) {x,y))n, 


(T( i ){x,y)^ n (x 1 ,y 1 )^- n {x 2 ,y 2 ))c 

(®n(xi,yi)®-n(x 2 ,y 2 ))c 


(15) 

and determines the weights of the 2d GCFT twist fields 
$n(x,y) as h L = H (1 - £) and h M = (l - ^)- 
Arguments analogous to the ones for 2d CFTs then per- 
mit us to infer 


Tr Pa = k n {®h L ,h M (xi,yi)$ hL , hM (x 2 ,y 2 ))c 


— k n X 12 1 


exp 


cm , 1 , y 12 

— (n- -) — 

L 12 n X12 J 


(16) 


for some constants fc n . We fix Aq conveniently. 

We derive EE as before from a limit of Renyi entropies, 
Sa = — lim„^! -J^Tr p^, and finally obtain 


£GCFT 2 = C Jl ln 

6 a 


CM 

~&T X 


(17) 


where t x = 0 x 12 , £ y = ai/ 12 , and a is again the lattice size 
that regularizes ultra-violet divergences. The EE formula 
(17) is one of our main results and completes the first 
goal, namely to calculate EE for a 2d GCFT. Note the 
similarity of the first term in the GCFT EE with the 
expression for the CFT EE (1). This is expected, since 
in theories with vanishing cm the non-trivial part of the 
GCA (6) reduces to a single copy of the Virasoro algebra 
[31, 53], 

Generalization to finite temperature or finite size. We 
generalize now our main result (17), starting with the 
consideration of finite temperature T = /3 -1 . We fol- 
low again the lead of Calabrese and Cardy [4]. For this 
purpose, we utilize the fact that Tr behaves like the 
2-point function of some primary fields in a GCFT. We 
use the transformation laws of the fields from the plane to 
the cylinder elucidated in [44] to find after some algebra: 



^Pl2COth(^12) (18) 


This is the EE for a mixed state in the 2d GCFT at finite 
temperature /3 _1 . 


The result (18) has the expected behaviour under two 
extreme limits. The first limit, £ 12 , P 12 <C /?, gives back 
the answer on the plane (17) with £ 12 —> t x l a an( l P 12 — > 
iy/a. The second limit, £ 12 , p\ 2 /3, yields 

S' = ^(ci£r 2 +c M P 12 ) + y In /? + •■ ■ (19) 

where the ellipsis denotes temperature-independent and 
subleading terms in /3. To leading order the EE (19) is ex- 
tensive, i.e. , linear in the separations, which is completely 
analog to what happens in 2d CFTs [4]. The subleading, 
logarithmic, term is universal, depending only on cl- 

To find the entropy of a subsystem of length (£ 12 , P 12 ) 
in a finite system of length L in its ground state, we use a 
similar analysis, replacing f3 — > L (and orienting the cuts 
along which the twist operators lie perpendicular to the 
previous case, which results in some changes). We then 
obtain a result for EE in systems of length L : 

S = y ln sin (^ 12 )) + cot (^£ 12 ) (20) 

Towards higher dimensions. One of the great advan- 
tages of GCFTs is that the infinite dimensional symme- 
try extends to all spacetime dimensions. These have been 
realized in physical systems e.g. hydrodynamics [28] and 
more recently in Galilean electromagnetic theories [55]. 
It is plausible that the constructions above would work 
similarly for higher dimensional cases. 

Holographic EE in GCFT. In the remainder of our 
work we focus on the second goal, namely to derive our 
results for EE holographically. As a first step, we summa- 
rize aspects of proposed 3d gravity duals to 2d GCFTs. 
A minimal requirement for a putative gravity dual is that 
its asymptotic symmetries match the symmetries of the 
corresponding QFT. Flat space holography in 3d real- 
izes this requirement [29, 30]. (See [56] for a different 
approach to holography and entanglement in flat space.) 

Flat space geodesics. Following the Ryu-Takayanagi 
prescription [8] one can determine EE holographically 
in 3d (Einstein) gravity by computing the length of 
geodesics on equal time slices. To determine the rele- 
vant geodesics we consider the most general zero-mode 
solutions of 3d flat space Einstein gravity, given by the 
line-element [34] 

ds 2 = M d u 2 — 2 du dr + J du d<p + r 2 dtp 2 (21) 

with mass parameter M and angular momentum param- 
eter J. For M = J = 0 the solution is known as null 
orbifold [57-60]. For M = — 1 and J = 0 the solution 
is (global) flat space. The more generic case M > 0 
and J/0 corresponds to flat space cosmology solutions 
[61, 62], which we will not consider in the current work. 

The central charges for flat space Einstein gravity are 
Cl = 0 and Cm = 3/Gjv, where Gn is the 3d Newton con- 
stant [30]. Therefore, the only term that remains in EE 
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for flat space Einstein gravity is proportional to i v [pi?\ 
on the plane at zero temperature [at finite temperature 
or for finite size], see (17) [(18) or (20)]. In the coordi- 
nates above this separation translates into the separation 
u\ 2 - This means that taking a constant u-slice we should 
obtain vanishing EE. Indeed, this is precisely what one 
obtains from calculating the length of such geodesics. 

While the conclusion above shows consistency of the 
holographic calculation with the GCFT results, it is more 
interesting to consider the general case where both cen- 
tral charges are non-zero. For this purpose, we no longer 
can consider geodesics of Einstein gravity solutions. 

Holographic EE from Chern-Simons formulation. 
Einstein gravity in 3d can conveniently be reformulated 
as a Chern-Simons gauge theory [63, 64] with the action 

I CS [A] = ^J(AAdA + ^AaAaA), (22) 

where k is the Chern-Simons level, A takes values in some 
gauge algebra g and (...) denotes a suitable invariant 
bilinear form on the gauge algebra. In flat space the 
gauge algebra is given by g = is 1(2), which is precisely 
the algebra (6) restricted to n = 0, ±1. The flat space 
connection describing spacetimes (21) reads [65] 

A = b~ 1 (d+a)b with b = e*"- 1 , (23) 

a = ^ Mi —M-i'j cki + (Li — — L_i — ^ ' 

Concurrent with the Grassmanian interpretation of [43] 
we refer to the generators L n ( M n ) as “even” (“odd”). 

In [66-68] it has been argued for AdS 3 that using the 
Chern -Simons formalism holographic EE of a given inter- 
val at the boundary of AdS is given by (the logarithm of) 
a Wilson line attached at the endpoints of this interval. 

S EE = ~\n{W n (C)} 

W n {C) = Tr TC ( Ve-Tc A ) = J VUe~ s{u ’ A)c 

Here 1Z denotes a suitable representation of the symme- 
try algebra at hand, V stands for path-ordering, C de- 
notes the contour of the Wilson line, and U refers to a 
topological probe exploring the bulk geometry. 

For flat space we assume a split of the probe U and the 
action S{U\A)c into even and odd parts labeled by L 
and M, respectively, i.e., U = U L U M , and S{U\A)c = 
Sl{Ul]Al)g + Sm(Um', Am)c , where A = Al + Am- 
Using this assumption the even action Sl is given by 

Sl = ^d s ((P i p i U L )f/^ 1 > i +A L ((P|) i -C 2 i )) (25) 

where Dl = 8 s +Al , Pl is the canonical momentum con- 
jugate to Ul, and A l is a Lagrange multiplier that sets 
the quadratic Casimir to the right value Clf — c|/288. 
Note that at this stage we do no longer assume that we 


are necessarily dealing with Einstein gravity, where cl 
vanishes; an example [31] for a theory with cl ^ 0 is 
topologically massive gravity [69] whose AdS EE was de- 
rived holographically in [68]. The brackets (■ ■ -)l denote 
a suitable invariant bilinear form satisfying (Pl, Pl)l = 

pL-ipLi + IpLopLo _ pLipL-i The odd action Sm 

(and all quantities therein) are obtained from the even 
action (25), replacing everywhere indices L — > M . 

We determine now the equations of motions from these 
two actions and find solutions to them, using the so- 
called “nothingness” trick: we first solve the equations 
of motion for the trivial case Al = Am = 0 and de- 
termine U^m and Plm- Then we use a (large) gauge 

transformation U(s ) = g(s)U^U^ g~ 1 (s), where the 
group element g(s) is in the component continuously 
connected with the identity to obtain solutions for non- 
trivial gauge fields Al,m [72], The actions Sl and Sm 
then reduce on-shell to = —2IS.oll,mC^ M with 

A a L ,M = a L ,M(sf) - ol,m( 0) determined by U^ M . 

The only thing left to do is to choose suitable boundary 
conditions for the topological probe U(s) at the begin- 
ning s = 0 and the end s = Sf of the boundary interval 
and to determine A oll.m- The reason why it is sufficient 
to determine A oll^m is that in the semiclassical limit 
/ VUe~ s{ - u ^ c ~ e -S(u,A)c and t h U s 

See = — 2A a L C% - 2A a M C^ . (26) 


We have found suitable boundary conditions, whose de- 
tailed discussion will be presented elsewhere [70]. 

From the non-trivial connection (23) and its split into 
even and odd parts, Al + Am = gdg~ x with g = 
b~ 1 e~f aidx , one can now readily obtain the EE for 
various interesting spacetimes. We collect these results 
below. We define i u = Uf — uo and ^ = ipf — ipo; 
the quantity a is again a regularizing cut-off. Applying 
the procedure above yields the EE of the null orbifold 
(M = J = 0), to be compared with (17). 


cl , Cm ^ u 

-6 ln ^ + ~6 'J, 


(27) 


Similarly, for global flat space (M = — 1, J = 0) we 
obtain the EE [to be compared with (20) for L = 27r] 

See = y ln ( 2sin y) + y^cot^A . (28) 

The holographic EE results above are in precise agree- 
ment with the ones obtained on the field theory side. 
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